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I. INTRODUCTION 



At present two deep inelastic scattering (DIS) sum rules, the Gross-Llewellyn Smith 
(GLS) |l[ and the polarized Bjorken sum rule 0, give the possibility of extracting the value 
of the strong coupling constant as from experimental data, in particular, at low momentum 
transfer, Q, down to 1 GeV (see, e.g., Refs. [0,3). Comparison of these values of as with 
other accurate as values, such as those obtained from the decay widths of the r-lepton and 
the Z-boson into hadrons, is an important test of the consistency of QCD. Clearly what are 
required are reliable theoretical relations between physically observable quantities and the 
strong coupling constant as- At low Q 2 scales, there are significant theoretical uncertainties, 
which come, firstly, from a truncation of the series obtained from perturbation theory (PT), 
which leads to a significant renormalization scheme (RS) ambiguity, and, secondly, from 
poorly understood nonperturbative effects (see, e.g., Refs. [|]|| for a review). In this paper 
we apply the method proposed in Refs. |7]|| (see also Refs. H[IIJ), the so-called analytic 
perturbation theory (APT), to study the GLS sum rule, continuing our investigation of the 
APT approach initiated in Refs. |TT|JT2]| . This method takes into account the fundamental 
principle of causality, which in the simplest cases is reflected in the form of Q 2 -analyticity 
of the Kallen-Lehmann type. The standard renormalization group resummation violates 
this required analytic structure, and unphysical singularities such as a ghost pole appear. 
The APT approach maintains the causal properties and removes all unphysical singularities 
by incorporating nonperturbative terms (for the number of flavors occurring in nature) 
to secure the required analytic properties. These nonperturbative terms can be presented 
as power corrections and their appearance is not inconsistent with the operator product 
expansion 

It is familiar that the QCD correction to the GLS sum rule is expressed as an expansion 
in powers of as that allows one, in principle, easily to obtain a value for the running cou- 
pling constant. Before going into detailed considerations, let us demonstrate the difference 
between the standard PT and APT running coupling constants in the two-loop level. The 
two-loop analytic running coupling constant can be written in the form of a sum of the 
standard perturbative part and additional terms which compensate for the contributions 
of unphysical singularities, the ghost pole and a cut arising from the "log-of-log" depen- 
dence. The difference may be transparently shown by using the approximate formulae for 
the two-loop analytic running coupling given in Ref. M 



«apt(Q 2 , A, n/) ~ a PT (Q 2 , A, 7i/) 
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(1) 



where (3q = 11 — 2n//3 is the one- loop coefficient of the /5-function corresponding to rif 
active quarks, and, for four active flavors, C\(nf — 4) = 0.0396. The PT running coupling 
constant is obtained by integration of the renormalization group equation with the two-loop 
/3-function. 

The difference between the APT and PT functions is illustrated in Fig. [I] over a wide 
range of Q 2 , 1 < Q 2 < 10 4 GeV 2 . The curves in the figure correspond to various values of 
the QCD scale parameter A for four active flavors. Fig. [I] shows that the expression (|l]), 
represented by dotted lines, approximates the exact two-loop analytic coupling (solid lines) 
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rather well for Q > 2A. Comparing the Q 2 evolution of the QCD running coupling constant 
obtained from APT to that given by PT (dashed lines), one can see that the difference 
between the shapes of the APT and PT running coupling constants becomes significant at 
low Q 2 -scales, Q 2 < lOGeV 2 . This fact stimulated applications of the modified perturba- 
tion theory with correct analytic properties, APT, for various physical processes (see, e.g., 
Refs. ||II1 , |14}1 ). Here, we consider the GLS sum rule in the framework of the APT approach. 
As has been demonstrated in Ref. JT5] by using the Deser-Gilbert-Sudarshan representation 
for the virtual forward Compton amplitude |16| (see also Ref. [|17|), the moments of the 
structure functions are analytic functions in the complex Q 2 -plane with a cut along the 
negative real axis. On the other hand, the conventional renormalization group resummation 
does not support these analytic properties and the influence of requiring these properties 
to hold in the DIS description has not been studied. Here, we perform this investigation, 
by applying the APT method, which gives the possibility of combining the renormalization 
group resummation with correct analytic properties of the QCD correction to the GLS sum 
rule. In Sec. II we start by describing the GLS sum rule in the PT and APT approaches 
and compare the Q 2 evolution of the APT and PT predictions. In Sec. Ill we consider in 
detail the RS dependence of the results. Summarizing comments are given in Sec. IV. 



II. THE GLS SUM RULE WITHIN APT 



The GLS sum rule predicts the value of the integral over all x of the non-singlet F3 
structure function measured in neutrino- and antineutrino-proton scattering 



— [ dx 
2 Jo 



F^(x,Q 2 ) + F^(x,Q 2 



(2) 



In the quark-parton level, which is appropriate for Q 2 — > 00, the GLS sum rule should equal 
three. Therefore, for fixed Q 2 the integral (0) can be conveniently written as 



£gls(<2 2 



1 - A GLS (Q 2 



(3) 



where the QCD correction, Aqls, m principle, contains perturbative and nonperturbative 
parts. To begin, we concentrate on the perturbative contribution to Agls ; considering in 
turn standard PT and APT methods and postponing until later a discussion of the possible 
influence of higher twist (HT) effects, which remain poorly understood. In this connection 
it is interesting to note from the result of the fit to the revised CCFR data |18| presented in 
Refs. [^] that the higher-twist contributions are small in the region Q 2 > 1 GeV 2 and can 
have a sign-alternating character. 

The standard perturbative part of the GLS sum rule correction is known up to the three- 
loop level for massless quarks in the MS-like renormalization schemes with the number of 
active quarks rif fixed, 



a pt 

GLS 



apT +d 1 (RS,n / ) (^)\d 2 (RS,n f ) ^ PT ^ 



7T 



V 7T 



7T 



(4) 



The PT running coupling constant «pt(Q 2 , RS, n/) is obtained by integration of the renor- 
malization group equation with the three-loop /^-function. The coefficients d\ and c?2 are 
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given in Ref. |2(J in the MS renormalization scheme. Therefore, the perturbative QCD cor- 
rection to the GLS sum rule is represented in the form of a power series in apx and, at 
first glance, the value of «pt can be easily extracted if the value of A GL s(Q 2 ) is experimen- 
tally known. In the region Q 2 > 1 GeV 2 , one believes PT with its renormalization-group 
improvement is still valid. We should note that different regions of the x integration in 
Eq. (0) at fixed values of Q 2 , in principle, correspond to different numbers of active quarks, 
nf, arguments have been given in Ref. PI] to select one or another value of n/. Experi- 
mental measurements of the GLS sum rule are made in a region of x where one believes 
that four light flavors are relevant. At present, there is no regular and consistent method of 
including threshold effects for the GLS sum rule. In the following analysis, we will first take 
rif = 4 to obtain results in the standard MS renormalization scheme, and then consider RS 
dependence. 

The QCD correction with correct analytic properties can be written in the form of a 
spectral representation 



Agls(Q 2 ) 



1 

7T JO 



da 



a 



Q- 



(5) 



where we have introduced the spectral function, which is defined as the discontinuity of 
A GLS (Q 2 ): f?(cr) = Disc {A G ls( — °" — /2i • If we calculate now the spectral function 
g(a) perturbatively, we get an expression for Agls(Q 2 )) which has the correct analytic 
properties and therefore no unphysical singularities. Consequently, we write the three-loop 
APT approximation to Agls(Q 2 ) as follows 



Agls T = s ( aMQ 2 ) + d x s^UQ 2 ) + d 2 *)SWQ 



r(i) 



(2) 



f(3) 



(6) 



where the coefficients d\ and d 2 are the same as in Eq. (f|) and the functions £apt(Q 2 ) are 
derived from the spectral representation and correspond to the discontinuity of the fc-th 
power of the PT running coupling constant 



<5apt(Q 2 
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(7) 



The function S^l T (Q 2 ) defines the APT running coupling constant, «apt(Q 2 ) = ^apt(Q 2 ) > 
which in the one-loop order is given by 



(i) 



«apt(<5" 



47T 
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In (Q 2 /A 2 



+ 



Q 2 /A 2 



(8) 



As can be seen from Eq. (|6|), the first term of the expansion is a apt /tt, but the following 
terms are not representable as powers of oapt unlike in the PT case. There are approximate 
expressions, like Eq. ([I]), for higher loop corrections in Eq. (^]), which have rather simple 
forms, which can be derived by using a method of subtracting unphysical singularities |f22| . 

(2) 

For instance, for 5Xpt? ^ ne approximate formula is 



S {2) (z) 
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Q 2 /A 2 



(9) 
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where C\ = 0.0273 for four active flavors and (3\ = 102 — 38n//3. 

To illustrate the difference between the convergence properties of the PT expan- 
sion (H) and the APT series (||) we use the recent result of the CCFR Collaboration 
(CCFR'97) mil: S GLS = 2.47 ± 0.09 at Q 2 = 3.1 GeV 2 , which is consistent with the 
result of a previous CCFR analysis (CCFR'93) J2§, S GLS = 2.50 ± 0.08 at Q 2 = 3 GeV 2 . 
The central value of Sqls corresponds to the value of the QCD correction Aq£s = 0.1767 
and the successive terms of the PT series respectively constitute 65.1%, 24.4% and 10.5% 
of the total. At the same time, the corresponding contributions to the APT series (^) make 
up 75.7%, 20.7% and 3.6% of the total. The convergence of the APT series seems to be 
somewhat better behaved than is that of the PT expansion at such small Q ~ 1.76 GeV . 

The same may be seen from Fig. |2|, where Sols is shown as a function of the QCD 
running coupling constant as in the PT and APT approaches. As outlined above, in the 
PT case, the function Sgls is an explicit function of the PT running coupling constant and 
in the one-loop approximation is represented by a straight line in Fig. ^, as a parabola in 
the two-loop case, and as a cubic curve in the three-loop one. At sufficiently large values of 
as ~ 0.4, the difference between the 1-, 2-, and 3-loop PT predictions becomes large. An 
inclusion of the higher-twist term with the value recommended by the Particle Data Group, 



A HT = (0.09 ± 0.045) IQ 2 |2§ (see Refs. |Z7|,f28J for additional details), also significantly 
changes the behavior, as is apparent from the figure. Note that the coincidence of the one- 
loop PT and APT curves in Fig. [| does not mean that the PT and APT approaches are 
physically identical, this is simply a matter of the linear form of the one-loop approximation; 
the behavior of the PT and APT running coupling constants are rather different [see Eq. (§)]. 
In the APT case, the contribution of the higher loop corrections is not so large as in the 
PT one and the corresponding curves in Fig. |2| are quite close to the linear function, and, 
especially, there is very little difference between the 2-loop and 3-loop results. The horizontal 
lines in Fig. ||] correspond to central values from experimental data at different low values 
of Q 2 : S'gls (1-7 GeV 2 ) = 2.13 ± 0.46 [F2TJ, Sgls(3 GeV 2 ) = 2.50 ± 0.08 @, 5 GLS (5GeV 2 ) = 
2.63 ± 0.09 |23|,|24|]. The intersection of one of these lines with a given theoretical curve 
gives the value of as for that value of Q 2 in that theoretical description. It will be seen 
that stability of the theoretical curve at Q 2 ~ a few GeV 2 is required in order to extract a 
reliable value of as or the QCD scale parameter A. 

Consider the Q 2 evolution of the GLS sum rule. In Ref. p0|, the low-Q 2 dependence of 



GLS sum rule has been evaluated by combining xF% measurements of the CCFR'93 data 
with data from other DIS experiments. Preliminary updated analysis with CCFR'97 data 
was presented in Refs. |p3j , p4}] . The analysis of the GLS sum rule based on the Jacobi poly- 
nomial expansion has been given in Ref. JJTJ, and for new CCFR'97 data has been examined 
The very recent CCFR/NuTEV result for the GLS integral is presented in 



f. 


3; 
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Refs. H||. The three values of the GLS integral at Q 2 = 1.26, 2.00, and 3.16 GeV 2 are 
in good agreement with the old result [[SIJ, but the value at Q 2 = 5 GeV 2 is larger, which 
corresponds to a smaller value of the QCD correction to GLS sum rule, although consistent 
within errors (see also Fig. 4). Note that the value SqlsW 2 = 5-01 GeV 2 ) = 2.776 |53J gives 
a very small value of the scale parameter Ap^ -4 '* = 89 MeV. The choice of normalization 
point influences the value of the parameter A, but does not change the general picture of 
there being a difference between the APT and PT results for low Q 2 scales. To study this 
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difference we take the relevant value of A^ 4 ) = 300 MeV. That practically means the nor- 
malization on the value S GLS (Q 2 = 3.16 GeV 2 ) = 2.547 that gives A^, = 303 MeV. 
This value of the scale parameter A is quite realistic and agrees well with the results of the 
fits to the structure function xF% from the CCFR'97 data, for example, with the two-loop 
result 337 ± 28 MeV pf and with the three-loop value 308 ± 34 MeV @. Since the differ- 
ence between the PT and APT forms of the QCD corrections is of order A 2 /Q 2 , both these 
functions will coincide in the asymptotic region, where the perturbative approximation is 
valid. 

The comparison between the Q 2 evolution in the APT and PT approaches is shown in 
Fig. |3], where the QCD correction to the GLS sum rule is plotted for the perturbative part 
(solid curve for APT and dashed for PT approach) and separately for the HT term given 
by two different estimations for the coefficients in the form (0.09 ± 0.045)/Q 2 (dash-dotted 
line) and (0.16 ± 0.01)/Q 2 (dotted line) taken from Refs. [P^PSfl , respectively. This figure 
demonstrates that there is an essential difference between PT and APT evolutions for low 
Q 2 . Instead of a rapidly changing function with unphysical singularities as occurs in the PT 
case, we get a slowly changing function in the APT approach. 

In Fig. H we plot the full contribution to the QCD correction Aqls with the perturbative 
part, calculated in the PT and APT approaches, and the HT part taken in the form as dis- 
cussed above f26fl . We also plot the data, indicated by squares, of the recent CCFR/NuTeV 
analysis [0,0 where the Q 2 behavior of the GLS integral has been evaluated at values of 
Q 2 between 1.26 GeV 2 and 12.59 GeV 2 , as well as older data. The dotted curve represents 
the PT result without HT effects. At high Q 2 scales, the PT and APT results agree closely 
with each other, both including the HT effects and without them, whereas at low Q 2 scales, 
the difference between the PT and APT behaviors becomes significant. At the same time, 
for Q 2 > 2 GeV 2 , the PT prediction without the HT term practically coincides with the 
APT curve including the HT term. Cancellation between the additional APT terms beyond 
the standard PT prediction and the HT terms can explain the fact that attempts to extract 
the HT effect from the CCFR'97 xF 3 experimental data give a rather small value, even 
poorly determining the sign of the HT contribution (see, for detail, Ref. |]19|). Note also 
that in the analysis of Refs. fl23"l , |2"4|| the HT contribution to the QCD correction for the GLS 
sum rule is taken to be given by the form (0.09 ± 0.05)/Q 2 ; however, in the subsequent 
papers p3|-|3~5|j , the HT term has been taken to be smaller, (0.05 ± 0.05)/Q 2 . Thus, the 
comparison of the APT and PT predictions with experimental data, as is demonstrated in 
Fig. [|, cannot give a definite conclusion since there are large experimental errors and the 
value of the HT corrections is very uncertain. 

We have considered the Q 2 evolution of the GLS sum rule in the customary MS renormal- 
ization scheme. In general, it is not sufficient to obtain a result in some scheme, but rather 
it is important to study its RS stability over some acceptable domain. The RS dependence 
of the GLS sum rule based on the PT approach has been studied in Ref. ||21|| . In the next 
section we consider the RS stability of the APT results. 
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III. RENORMALIZATION SCHEME DEPENDENCE 



A truncation of a perturbative expansion leads to uncertainties in the theoretical predic- 
tions arising from the RS dependence of the partial sum of the series. At low momentum 
scales these uncertainties may become very large (see, for example, an analysis in Ref. |p8|| ). 
A physical quantity, in our case the QCD correction to the GLS sum rule, has to be invariant 
under a change of RS, when the coupling constant transforms as follows (a = as/V) 

a' = a(l + Via + v 2 a 2 H ). (10) 

In the new RS, the QCD correction Aqls is represented as 

Agls = a' (1 + d[a' + d' 2 a' 2 ) , (11) 

where the coefficients d\ and d 2 are RS dependent and have the transformation law 

d' 1 = di-v ll (12) 
d' 2 = d 2 — 2(di — vi)vi — v 2 . 

In the three-loop level there are two process-dependent invariants under the RS-transfor- 
mation (p~0|) [p^j 

Pl = 2 A2 X ' ( ' 

p 2 = b 2 + d 2 - b x di - d\ , (14) 

where b = fio/2, b\ and b 2 are the coefficients of the three-loop /3-function 

(3{a) = fi 2 -^ = -f3 a 2 {l + ha + ha 2 ) . (15) 
o /i z 

The coefficient b 2 , as well di and d 2 , depends on RS, and, because of the first equation 
in (|T2|), the scale parameter A transforms as follows [ 4T)| : A' = A exp(vi/b). In terms of 
the scale parameter A^, the three-loop perturbative running coupling in any RS obeys the 
equation 

h In ( ^- ) =df E -d 1 + $(a, b 2 ) , (16) 



2 ~ I A 2 — 

MS 



where 



$(a,6 2 ) = --6 1 lni±^ (17) 
a po a 



+ b 2 



dx 



l + 6ix)(l + 6 1 x + 6 2 a; 2 ) 



Thus, any RS taken from the MS-like schemes can be characterized by two parameters, 
which we choose here to be di and b 2 . To calculate the second RS invariant one can use the 
coefficients in the MS scheme, which for nj = 4 are di = 3.25, b 2 = 3.05 . 
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In the framework of the conventional approach there is no solution of the RS dependence 
problem apart from calculating many further terms in the asymptotic PT expansion, and 
there is no fundamental principle upon which one can choose one or another preferable 
RS. Usually, one uses a class of 'natural' or 'well-behaved' RSs, which are defined by the 



so-called cancellation index criterion [|3]], according to which the degree of cancellation 
between the different terms in the second RS-invariant p 2 are not too large, as measured by 
the cancellation index 

C= T ^ 1 (\b 2 \ + \d 2 \ + d 2 + \d 1 \b 1 ) . (18) 

\P2\ 

By taking some maximum value of the cancellation index C max one can investigate the 
stability of predictions for those RSs with C < C max . In the case of the MS-scheme the 
value of the cancellation index for the GLS sum rule is C-^g = 95 . Taking into account 
that the MS scheme is commonly used, we will consider this value as a boundary for the 
class of 'natural' schemes. In addition we will compare our results with predictions given by 
optimized schemes based on the principle of minimal sensitivity (PMS) |39| and the method 
of effective charge (ECH) Q. (See also applications in plp3|, p| ). 



It should be stressed that the large value of the cancellation index C^jg for four active 
quark flavors does not mean that the coefficients from which the second RS invariant is 
constructed are huge, only that p 2 is small, p 2 = —0.32 . In this special case the value 
of p 2 is very sensitive to the value of n/; for example, for n/ = 3, this invariant becomes 
p 2 ni = 4.2 and the cancellation index becomes significantly reduced. At the same time, 
the absolute values of the coefficients that construct the second RS invariant (III) are smaller 



for rif = 4 than for nj = 3, however, C^~^ — 10 C^~ 3 \ Perhaps, in this situation, it 
is more convenient to introduce as the corresponding index, which will define the class of 
'natural' RSs, the sum of the absolute values of the terms in p 2 , without the denominator in 
Eq. (TBp. The numerical parameter introduced in such a way will practically define the same 
region in the (d\, b 2 )-pla.ne, without changing significantly with changing rif, and, therefore, 
more adequately describing the situation. 

Let us briefly discuss results of the ECH and PMS approaches. For the ECH scheme 
the parameters are <if CH = c?f CH = and, therefore, 6f CH = p 2 . The transformation 
from the MS scheme to the ECH scheme is performed with the parameters in Eq. ([TDD 
v ech = d Ms = 3 25 anc i v ech = d Ws = u 20 . The system of equations for getting the PMS 
optimal prescription consists of four equations. These are Eqs. (13) and (|16|) complemented 
by the following two equations: 

3d 2 + 2d x b x + b 2 + (2d l b 2 + 3d 2 &i) a + 3d 2 b 2 a 2 = , 

(1 + b x a + b 2 a 2 ) (1 + 2d x a + 3d 2 a 2 )I(a, b 2 ) - a = . (19) 

where 

IM= l + + <20) 

In the PMS optimization procedure the coefficients d\, d 2 and b 2 become Q 2 dependent and, 
therefore, have to be adjusted for each different value of Q 2 . 
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We present our results in Fig. |5[ where the QCD correction Aqls is plotted as a function 
of Q 2 in different RS. The schemes A and B have the same cancellation index as in the MS 
scheme, C = 95, and are defined by the following values of the parameters: df = —3.93, 

= and = —1.0, &f = 13.7. The dotted curves present the ECH and PMS predictions, 
which are very close to each other. This is because the PMS coefficients <if MS (Q 2 ) and 
<i PMS (Q 2 ) are found for four active quark flavors to be small numerically in the region under 
consideration, Q 2 > lGeV 2 . Behaviors of these coefficients as functions of Q 2 are shown 
separately in Fig. ||. 

As it has been mentioned above, there are no strong arguments to fix the definite number 
of active quark flavors for the integral (0) to be rif = 4. Therefore, we have also considered 
the case n/ = 3. The ECH and PMS prescriptions gave again results that are very close 
to each other. However, the reason is no longer that the PMS coefficients d PMS (Q 2 ) and 
g?2 MS (Q 2 ) are very small, but that there is a strong cancellation between the last two terms 
in the expression for the QCD correction A GLS = a PMS (l + <5 PMS ) with 5 PMS = d PMS a PMS + 

a PMS j This fact is demonstrated in Fig. [7], where (a) shows the quadratic (solid line) 
and cubic (dotted line) contributions to A G ls for = 3. The sum of these contributions is 
shown in Fig. |7J (b) for nf = 3 (solid line) and for nj = 4 (dashed line). 

The APT predictions for the MS, A, and B schemes turn out to be stable for the whole 
interval of momenta and appears as a wide solid line in Fig. |5|. It should be stressed that 
the optimized RS constructions in the framework of PT do not maintain the correct causal 
properties, however, to avoid this difficulty, they can be modified in the sense of APT. In 
Fig. |5] we plot the analytic ECH result, which, of course, is different from the perturbative 
one and practically coincides with the other APT curves for the MS, A, and B schemes. 

The sensitivity of the PT predictions to the RS dependence can be reduced by applying 
the Pade approximation (PA) |35|] . However, in the case under consideration, this approach 
can lead to some difficulties. In Fig. [| we plot two possible Pade approximants [1/1] and 
[0/2] in the three-loop order for Aqls versus as- This figure demonstrates that the PA leads 
to unphysical singularities in the region of small momenta that violate the causal properties 
and prevents us from considering the method of PA as a systematic approach. In contrast, 
the APT results for the MS, A, B, and ECH schemes lead to stable predictions for the whole 
interval of momentum, which practically coincide with each other, and, therefore, the APT 
description turns out to be practically RS independent. 



IV. SUMMARY AND CONCLUSION 

We have considered the GLS sum rule by using the APT approach, which maintains 
the causal properties and, in contrast to PT, does not lead to any unphysical singularities. 
Taking into account the analytic properties of the moments of the deep inelastic structure 
functions, we have obtained an analytically improved theoretical description of the GLS 
sum rule. We have shown that the convergence properties of the APT expansion are better 
than are those in the standard PT description. We have demonstrated further that there 
is a significant difference between the Q 2 evolution in the PT and APT approaches for low 
Q 2 scales. We have also analyzed the GLS sum rule taking into account the higher-twist 
contribution. The APT power corrections have a sign opposite to that of the typically-used 
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higher-twist term p6] and, numerically, there is an effective cancellation between these two 
corrections. 

In this paper we have found that in the framework of APT the theoretical prediction to 
the GLS sum rule becomes practically renormalization scheme independent starting from 
the three-loop level. A similar statement has been made for the e + e~ annihilation into 
hadrons |P| and for the semileptonic inclusive decay of the r lepton [47] . A corresponding 
analysis of the Bjorken sum rule is presented in Ref. JTJ]]. Thus, the APT approach gives 
a systematic method of reducing the RS ambiguity significantly and leads to practically 
unique predictions for physical quantities. It should be stressed that our analysis shows that 
there is serious doubt concerning the conjecture that, in spite of the large values of as at 
low Q 2 , the conventional three-loop PT predictions for the GLS sum rule are reliable (see, 
e.g., [0). The proximity of PT predictions in the MS scheme to predictions obtained in 
the optimal schemes seems to possess no significance for any process at small momentum 
transfer and cannot be considered a guarantee of RS independence of the results. 

Unfortunately, at present, the experimental situation for the GLS sum rule, with its 
large errors, does not allow us to come to a reliable conclusion that the APT description 
is preferable to the PT one. However, from the theoretical point of view, the remarkable 
properties of the APT approach, for example, the maintenance of causality and the higher- 
loop and renormalization-scheme stability for the whole interval of momenta, create a basis 
for preferring the application of this new technique. 
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FIG. 1. Comparison of Q 2 evolution of the two-loop running coupling constant in the APT and 
the standard PT approaches. 
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FIG. 2. S'qls w hh 1-, 2-, and 3-loop QCD corrections vs. the coupling constant. 
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FIG. 3. The theoretical predictions for Q 2 evolution of QCD correction to the GLS sum rule 
given by Eqs. (Q) and (|6|). The dash-dotted and dotted lines represent the HT corrections from 
Refs. " 



27,28], respectively. 
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FIG. 4. The APT and PT predictions for the QCD correction to the GLS sum rule together 
with experimental results. The Serpukov data |29| is denoted by a circle, the CHARM [36| data 
by a downward pointing triangle, the CCFR data from Ref. 37] by an upward pointing triangle, 
and the recent CCFR/NuTeV results from Ref. [34| by squares. 
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